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Abstract. We present a dynamical model of heavy quark evolution in the quark- 
gluon plasma (QGP) based on the Fokker-Planck equation. We then apply this model 
to the case of central ultra-relativistic nucleus- nucleus collisions performed at RHIC 
and estimate the component of J /ip production (integrated and differential) stemming 
from c-c pairs that are initially uncorrelated. 



1. Introduction 

According to the first experimental results on open charm production in nucleon-nucleon 
reactions at RHIC {^/s = 200 GeV), up to 40 c-c pairs could be formed in a central 
Au-Au collision at this energy. As noticed in this large number could result in the 
formation of additional charmonia as compared to a simple scaling of nucleon-nucleon 
prompt production||, due toc-|-c— + reactions happening later in the QGP. 
One would then observe a charmonia enhancement instead of suppression. However 
(i) D + I) ip + X reactions happening during the ensuing hadronic phase can also 
contribute to charm enhancement and (ii) such an enhancement does not seem to be 
observed experimentally (preliminary J/ ip data appear to be compatible with the usual 
suppression). All of this indicates that the models should be refined. In fact, the actual 
distribution of c and c in phase-space could have a large impact on these recombination 
processes, but is often not treated with the care it deserves (or is simply beyond the 
scope of some models). Here, we will complement the model introduced in [1 by a 
dynamical treatment of the heavy-quarks evolution in the QGP based on a Fokker- 
Planck equation. We also present and discuss the differential spectra of c quarks and 
charmonia at the end of the QGP phase. 

§ To whom correspondence should be addressed (pol.gossiaux@subatech.in2p3.fr) 

II For a more detailed discussion of these ideas, see the contribution of R.L. Thews to this conference. 
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2. The model 

Following |2lini, the c and c quarks distributions in the QGP are assumed to follow a 
Fokker-Planck (FP) equation in momentum space. 



df{f>,t) d 



A{p}f{p, t) + ^ {Bij{p)f{p, t)) 



(1) 



dt dpi 

The main justification for this hypothesis is that the heavy mass of these quarks implies 
large relaxation times as compared to the typical time of individual c + g c' + g' 
and c + q — > c' + q' collisions, whatever the momentum of the heavy quark. The drag 
(A) and diffusion (B) coefficients were evaluated according to jS El resorting to 
a power expansion "a la Kramers - Moyal" (KM) of the Boltzmann kernel of 2 ^ 2 
collisions. As realized in 0, the asymptotic distribution coming out of the FP evolution 
deviates from a Boltzmannian. This results from the truncation of the KM series. 
We have therefore slightly corrected the B coefficients in order to guarantee a correct 
Boltzmannian asymptotic distribution. We have also checked on various examples that 
this modification has little effect at small evolution time. This constitutes our reference 
set of FP coefficients. In order to circumvent our lack of knowledge on the c-QGP 
interaction (radiative processes, non perturbative aspects, etc.^), we will also consider 
alternative sets, obtained by multiplying the reference one by a numerical factor K. It 
is our hope that the gross experimental results will permit to fix an approximate value 
of K and that our model could then be used to predict finer aspects. In this respect, 
we view our model as a semi-predictive effective theory that could be useful to match 
the gap between the fundamental underlying theory (QCD) and experimental results. 

The FP coefficients depend on position and time (only) through the local 
temperature and velocity of the surrounding medium, that is assumed to be locally 
thermalized and described via hydrodynamic evolution. Therefore, we do need to 
evaluate explicitly all microscopic c + g — > c' + g' and c + q c' + q' processes when 
performing our simulations: only c, c and J/ip d.o.f. are considered. 

For the sake of simplicity, we have neglected all initial state interactions in this 
exploratory work. We have also allowed ourselves to take a simple factorized form of 
the initial c — c phase-distribution, i.e. /m(^c5 Pc! ^e, Pe) oc TA{rc,±)TB{rc,±)5^^\rc — 

rc)S{Zc) X /in(Pc) X finiPc) with /in(Pc) OC /in(Pc,±) X /in(l/c), whcrC /in(?/c) haS bcCU 

chosen according to and /in(Pc,±) according to D transverse- momentum spectra at 
mid-rapidity in nucleon-nucleon reactions. 

J/i/j mesons will be assumed to exist up to a dissociation temperature T^iss that 
we let as a free parameter. Under this temperature, they behave as if they were in 
vacuum. According to recent lattice calculation, Tdiss could be substantially above the 
transition temperature. Apart from the prompt component, J/ip are assumed to be 
formed and destroyed in the QGP via the c-|-c ^ J/ip + g processes P]. In our 
simulation, these radiative dissociations are implemented via local dissociation rates 
that efficiently incorporate the microscopic processes. 
% Ultimately, transport coefficients should be evaluated using lattice calculations. 
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3. Results and discussion 



All results are for central Au+Au reactions at ^/s = 200 GeV. The initial number Nc 
of c and c quarks is taken to be 40, with ^^''jfy~^^ ~ 9. Unless specified differently, the 
temperature and velocity profiles have been evaluated with the hydrodynamic model of 
Kolb and Heinz jTj . In order to address most clearly the effect of the uncorrelated J/ ip 
production, the prompt J/ip component has been set to zero. 




Figure 1. differential productions of c quarks at the end of the QGP phase: On 
the left panel, the rapidity distribution is presented relative to the one at initial time. 
On the right panel, transverse momentum spectra are displayed for Kolb and Heinz 
hydrodynamic evolution of QGP (full lines), as well as for some Bjorken hydrodynamic 
evolution with no transverse expansion (broken lines). 

We show the differential productions profile of c quarks at the end of the QGP 
phase in figure Q for different values of K. Let us remind that typical relaxation times 
are oc K^^. Thus, K = corresponds to the free streaming case while K —* +oo 
corresponds to instantaneous equilibration of heavy quarks with the QGP. 

On a broad window around mid-rapidity, the final profile deviates merely from the 
initial one. One just notices a small broadening for large values of K. Indeed, heavy 
quarks locate themselves very early in a slot of matter going at the same rapidity. The 
drag force thus vanishes and the heavy quarks conserve their rapidity for the rest of 
the evolution. However, diffusion processes can still happen and bring heavy quarks 
towards the tail of the rapidity distribution, leading to some rapidity broadening. 

As for the transverse momentum distributions, the noticeable fact is the "cooling" 
of c quarks transverse motion by the QGP. One observes significant deviations from 
the free streaming case for > 5 only. For large values of K, distributions tend to a 
Boltzmannian at small and intermediate values of p±. As the distributions associated 
to different values of K interpolate smoothly between those two limiting cases, it might 
be possible - once initial and final state effects are included - to "read" the value of K 
from experimental spectra of D mesons. 

We present the differential production of J/ip at mid-rapidity on the left panel of 
figure El Tdiss appears to have a much larger influence on this quantity than K. Once 
the initial number of c — c is fixed more accurately, it might be possible to assess the 
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value of Tdiss and perform comparison with theoretical estimates. The right panel shows 
some typical predictions for J / ip transverse spectrum at mid-rapidity. In all cases - but 
especially for large K - the prompt component is harder that the uncorrelated one. 




Figure 2. Left: J/^ production at mid-rapidity at the end of the QGP phase with 
(full) and without (broken) radial expansion. Right: transverse momentum spectrum 
at mid-rapidity; plain lines correspond to A' = (hardest spectrum at large p_l), 1, 
5, 10, 20 & 40. The chain curve represents the best fit of the K = AO spectrum by a 
Boltzmannian, while the broken curve represents the prompt J /^p component. 



4. Conclusion and outlook 



We have presented a model that copes efficiently with dynamical evolution of heavy 
quarks and quarkonia in QGP. We have evaluated J/ip production at mid-rapidity 
as well as some differential spectra for different equilibration "strengths" K and for 
different dissociation temperatures. Hopefully, the comparison of such results with 
their experimental equivalent will ultimately allow to "measure" these parameters 
characteristic of the QGP phase. 
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